Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 6/22/04 



O 
O 
(N 



> 

O 

o 
o 

:^: 

Oh' 

6 



X 



SPATIALLY RESOLVED ULTRAVIOLET, Ha, INFRARED, AND RADIO STAR FORMATION IN M81 

K. D. Gordon\ p. G. Perez-Gonzalez\ K. A. Misselt\ E. J. Murphy^, G. J. Bendo\ F. Walter''*, 

M. D. Thornley^'', R. C. Kennicutt, Jr.\ G. H. Rieke\ C. W. Engelbracht\ J.-D. T. Smith^ 

A. Alonso-Herrero\ P. N. Appleton'*, D. Calzetti^ D. A. Dale**, B. T. Draine^, D. T. Prayer*, G. Helou**, 

J. L. HiNz\ D. C. HINES^"^ D. M. Kelly\ J. E. Morrison\ J. Muzerolle\ M. W. Regan^, J. A. Stansberry^ 

S. R. Stolovy*, L. J. Storrie-Lombardi**, K. Y. L. Su\ and E. T. Young^ 

Draft version February 2, 2008 

ABSTRACT 

We present Multiband Imaging Photometer for Spitzer (MIPS) observatfons of M81 at 24, 70, and 
160 fim. The grand design nature of M81 is clearly seen, showing two well resolved spiral arms 
containing numerous bright star forming regions. The MIPS images reveal a significant amount of 
cold dust associated with the spiral arms. We investigate the variation of the ultraviolet (UV), Ha, 
and infrared (IR) luminosities and star formation rate (SFR) indicators across the face of M81 using 
the MIPS images and archival UV and Ha images. For regions in M81, we find that UV and Ha SFRs 
(uncorrected for dust attenuation) are always lower than the IR SFR. The cause of this behavior is 
dust attenuation and/or using SFR calibrations appropriate for entire galaxies, not regions in galaxies. 
The characteristics of the dust attenuation for the regions indicate the dust grains and/or geometry 
are different from those in starburst galaxies. The behavior of the infrared-radio correlation in M81 
is seen to vary from the global average, with variations correlated with the morphology of M81. 
Subject headings: galaxies: individual (M81) — galaxies: spiral — galaxies: ISM — dust, extinction 



1. INTRODUCTION 

The galaxy M81 (NGC 3031) is one of the largest 
grand design spiral galaxies in the sky (14' x 27'). 
It is a member of the M81 galaxy group whose 
principle members are M81, M82, NGC 3077, and 
N G C 2976 and include a number of dwarf galax- 
ies (lAppleton. Davies. fc StephensonI 119811 iBovce et al.l 
120011) . Prominent tidal HI tails around M81 s how that it 
is int eracting with M82 and NGC3077 (Ym i. Ho. fc Ld 
119941). The Ceoheid distance to M81 is 3.6 Mpc 
IjFreedman et alJl200H) . The proximity of M81 has made 
it a favorite target for many investiga tions of galaxy 
jroperties from the X-ray to radio (e.g. Kaufman ct al. 
Devereux. Jacobv. fc CiarduUo 1995; Hill et al.. 



iTgtaiKong et al.l l2000:"S^rtz et al. 2003|. 

Many, well-resolved images of M81 from ultravio- 
let (UV) to radio wavelengths have been taken, with 
the notable exception of the infrared (IR). Previously, 
well resolved far-IR images of galaxies have only been 
possible for Lo cal Group g alaxies such as the Mag - 
ellani c Clouds llBraun et alJ Il998l iWilke et al.lj2003[) . 
M31 IjHaas et al.ll998|) . and M33 pippelein et al.l200^ . 
With the successful launch of the Spitzer Space Telescope 
(^Vcrncr ct al. 2004), it is now possible to map many 
large galaxies in the far-IR with good spatial resolution, 
good sensitivity, and in a reaso nable amount of time. 
Infrared Array Camera (IRAC, iFazio et all I2004J) and 
the Multiband Imaging Photometer for Spitzer (MIPS, 
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iRieke et al.l 12004) images of M81 were taken as part of 
the commissioning of Spitzer and showcased in the first 
Spitzer press re lease. The MIPS ima ges are presented in 
this paper and iWillner et all l)2004j) presents the IRAC 
images. M81 is one of the key galaxies in the Spit zer In- 
frared Nearby Galaxies Survey(SINGS, Kennicutt et alJ 
l2?ffll . 

We have chosen to investigate the variation of star for- 
mation rate (SFR) indicators and the infrared-radio cor- 
relation across the disk of M81 to highlight two of the 
questions which can be probed with these new observa- 
tions. We compare the MIPS observations with existing 
UV and Ha images to probe the behaviors of the IR, 
Ha, and UV SFR indicators across M81. Such compar- 
isons of SFR indicators have been made for global galaxy 
fluxes (eg.. Bell fc Kennicutt 2001; Kcwlcv ct al. 20o2 
lR,osa-Gonzalez. Terlevich. fc Terlevich 2002), but rarely 
has it been possible to resolve all three SFR indicators 
in a single galaxy. By studying the resolved behavior of 
these SFR indicators, a better understanding of their sen- 
sitivities to dust and stellar age can be determined. This 
will improve the accuracy of these SFR indicators for 
both resolved and global galaxy measurements. Radio 
emission is another pos sible SFR indicator, but it is no t 
as mature as the others |Haarsma et al."2000'; 'Bcll"2003') . 
As such, we have chosen to explore the more established 
infrared-radio correlation j Hclou et al. 1985; Condoa 
ITggir^amReddv. fc Condon..200L .Pierini et al...2003,) . 

2. DATA 

2.1. MIPS Observations 

Images of M81 at 24, 70, and 160 u m were obtained in 
a manner similar to that described in 'Engelb racht et alJ 
(|2004) for NGC 55. The MIPS images were reduced 
using the M IPS Instrument Team D ata Analysis Tool 
llGordon et al.. .20041 as described in lEngelbracht et all 
(120041). The uncertainties on the final absolute calibra- 
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Fig. 1. — The MIPS images of M81 are shown at full instrumental resolution. These resolutions are 6", 18", and 40" for 24, 70, and 
160 fiva, respectively. The field-of-view of the images is 20' X 30' which is more than large enough to cover M81. These mosaics have a 
total exposure time of approximately 80, 40, and 8 seconds per point at 24, 70, and 160 fim, respectively. 
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Fig. 2. — The MIPS and ancillary images of M81 are sliown convolved to have the 160 /im PSF which introduces a prominent Airy ring 
around bright point sources. The field-of-view of the images is 1 5' X 23'. Bright, foreground stars were masked in the UV and H« images 
prior to convolution. The ellipse denotes the region discussed in ilij.ll 
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tions are estimated at 10%, 20%, and 20% for the 24, 
70, and 160 /xm data, respectively. The final MIPS mo- 
saics are displayed in Fig. ^ The 70 and 160 fxm images 
exhibit linear streaks along the scan direction (roughly 
north-south) which are residual instrumental artifacts 
due to the time dependent responsivity of the Ge de- 
tectors. 

2.2. Ancillary Data 

Ultraviolet images of M81 were taken with the Ultra- 
yiolet Imaging Teles cope in far-UV and near-UV bands 
IjStecher et al.lll997tlMarcum et al.l l200Tl. The deepest, 
geometrically corrected, flux calibrated images for each 
band were retrieved from the Multimission Archive at 
Space Telescope. The Ha image of M81 was taken from 
iGreenawalt et a l. (1998). The radio image at 20 cm was 
created from archival Very Large Array (VLA) C-array 
data, using standard AIPS reduction techniques. 

2.3. Resolution Matching 

Each image of M81 was cropped to a common field-of- 
view (20' X 30') and convolved to match the point-spread- 
function (PSF) of the MIPS 160 /im image to allow for a 
consistent comparison of the images. We used a Fourier 
Transform technique for the convolution as described in 
lEngclbracht ct al. ( 2004) . All of the convolved images 
are shown in Fig. [51 The 160 /im PSF clearly has an 
Airy ring resulting in a prominent ring surrounding point 
sources in some of the images. 

3. RESULTS 

3.1. Morphology 

At 160 ^m resolution, the grand design nature of M81 
is very prominent in all but the R band image. The 
R band image displays a morphology dominated by the 
bulge, with only faint spiral arms present. At the other 
wavelengths, the spiral arms show a continuum of mor- 
phologies from being a string of point sources (Ha and 
24 /im) to exhibiting both point sources and an under- 
lying smoother distribution (UV and 70 /im) to fairly 
smooth arms with faint point sources (160 fiin and ra- 
dio). This difference in the appearance is related to the 
age of stars probed and the importance of dust attenu- 
ation. For example. Ha probes younger star formation 
than the UV images and, as a result, the UV images are 
smoother than the Ha image. Evidence that dust at- 
tenuation is affecting the arm morphology is seen in the 
southeast of M81 where a prominent piece of the galaxy 
is missing from the UV and Ha images (see ellipse in 
Fig. 12), but not from the IR or radio images. 

3.2. Star Formation Rate Indicators 

We investigated the behavior of the popular UV, Ha, 
and IR SFRs in M81. Ha and UV SFR indicators suf- 
fer from dust attenuation and probe stellar ages of 10 
and 100 Myr, respectively. The IR SFR indicator does 
not have an age bias, but does assume that all of the 
starlight from the star formation is absorbed by dust 
and re-emitted in the IR. As we use resolved images of 
M81 for all three SFR indicators, we can probe the bi- 
ases affecting each method in detail. We have extracted 
fluxes for specific regions in M81 from the 160 /im resolu- 
tion images, correcting the UV and Ha fluxes for a fore- 
ground Milky Way extinction with an E{B — V) ~ 0.08 




Fig. 3. — The regions for which fluxes were extracted are shown 
superimposed on the M81 24 fim image. The size of the extraction 
aperture is reflected in the size of the circle for each region. The 
average diameter of the apertures is 1.2 kpc. The color coding of 
the regions is given in Fig. |^ The green box in the lower right 
corner gives the sky aperture. At 70 fj,m, the sky noise contribu- 
tion to the uncertainty in the average aperture is ~2 mjy. This 
shows that the linear streaks are not affecting the 70 fim aperture 
measurements as all the fluxes are above 80 mJy. The field-of-view 
of the image is 20' X 30'. 



IjSchlegel. Finkbeiner. fc Davislll998^ . The regions are 
shown in Fig. |2| 

In order to convert measured UV, Ha, and IR lumi- 
nosities to SFRs, calibrations are needed. One way to 
determine such calibrations is to use stellar evolutionar; 
synthesis models (eg ., Fioc fc Roc ca-Volmerangc 199 
iLeitherer erallll999|i . iKennicuttI ((1998) has provided 
such calibrations with stellar age as sumptions appro- 
priate for entire galaxies. While the iKennicufrg 1)19981) 
calibrations will be invalid at some level for regions of 
galaxies, we have chosen to use them as they allow a 
physical basis for interpretation of luminosities. The to- 
tal IR flux was calculated using the TIR equation of 
iDale fc: Helo u ( 2002!) formulated specifically to determine 
the 3-1100 fim luminosity given MIPS measurements. 
All of the SFRs in this paper are uncorrected for dust at- 
tenuation as proper treatment of dust attenua tion is be- 
yond the scope of this paper feg.. lWitt & Gordon.2000.1 . 

The SFRs for the whole galaxy are estimated to be 
0.31, 0.38, and 0.89 Mq/vt from the integrated UV, Ha, 
and IR luminosities, respectively. The luminosities and 
SFRs determined in the regions defined in Fig.^lare plot- 
ted in Fig. ^ The arm regions display the largest area- 
normalized IR SFRs which is evidence that they have the 
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Fig. 4. — The luminosities for UV and Ha, ratiocd to the IR luminosity, are plotted versus the IR luminosity. The IR luminosity is 
normalized to the area in kpc^ of each extraction aperture. The corresponding SFRs and SFR ratios are indicated on the opposite axes. 
Representative uncertainties are shown. 
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Fig. 5. — The IR/1500 luminosity ratio is plotted versus the 
UV color (iJ2500 7^1500 ratio). The symbols are the same as in 
Fig. |3] The solid curve gives the linear fit [Lxir/\Li$oq = 
-2. 42+3. 81(L2500 7^1500 )] to the M81 data. The star symbols give 
the observati ons for a sample of starburst galaxies and the dotted 
line gives the tMeur.er.J^ckmaii. fc Calzetti 1 1999) relationship de- 
termined for starburst galaxies. Representative uncertainties are 
shown. 



highest SFRs. The ratio of Ha/IR luminosities is higher 
for the arm regions than any of the other regions. For 
aU the regions, the UV and Ha SFRs are always smaller 
than the IR SFRs and there is a weak trend towards 
smaller UV/TIR and Ha/TIR SFR ratios as the TIR 
SFR increases. There are three possible causes of these 
effects: 1) Dust is attenuating the UV and Ha measure- 
ments. 2) A fraction of the IR flux may not be associ- 
ated with recent star formation (i.e., dust he ated by old 
stars) . 3) The luminosity-SFR calibrations of lKennicutti 
l)1998t) are not appropriate for individual regions in a 
galaxy, where there may be transfer of energy between 
different regions (radiative transfer effe cts) and7or stellar 
age differences from those assumed bv iKennicuttI ()1998l) 
which will affect SFR indicators at different levels. We 
will investigate the relative importance of these factors 
in future papers but briefly discuss the first two points 
below. 

The association of significant IR emission with bright 
UV and Ha regions is evidence for significant dust at- 



tenuation. The dust attenuation is probed in Fig. [3 
where the IR/UV ratio and UV color attenuation mea- 
sures are plotted. The IR/UV ratio has been shown to 
be a good dust attenuation measure for galaxies and the 
UV attenuations for the extracted regions vary between 
0.5 and 3.0 mag for log (LrfB.7ALiRno) between 0.5 and 
1.5 l|Gordon et alJl200n(). The L250o/-^i500 ratio is also a 
dust attenuation measure, but as it is a UV color, i t is sig- 
nificantly more uncertain l)Witt fc GordonI I2000D . This 
flux ratio is very si milar to the more familiar dust at- 
tenuation measure IjMeurer. Heckman. fc Calzetti 1999). 
The starburst calibration of these measures is shown in 
Fig. [S] and does not describe the measurements of M81. 
De viations from the starburst calibration were also found 
by'B_ell et al." (2002") for Large Magellanic Cloud HII re- 
gions and by Kong et al. (2004) for global measurements 
of galaxies more quiescent than starbursts. This implies 
that the dust geometry and/or dust grain properties are 
different between starburst galaxies and star forming re- 
gions in normal galaxies. 

In studies of IR SFRs, it is common to assume that the 
cold dust is mainly heated by old stars, and thus the cold 
dust emission (160 /im) should not be used in determin- 
ing IR SFRs. The M81 images presented in this paper al- 
low us to probe this assumption for the M81 disk but not 
the AGN dominated nucleus. The morphology of the spi- 
ral arms progresses from point-like to fairly smooth from 
24 to 160 ^m. The Ha and UV bright regions are seen 
at 160 /im implying that young stars heat a fraction of 
the 160 /im emission. The morphology of the IR images 
best matches the Ha and UV images as opposed to the 
R band image (an old star tracer in M81), but this cor- 
relation is not enough to identify the dominant 160 fiui 
heating source, as the spiral density wave concentrates 
the dust in the arms. The smoother morphology of the 
UV images compared to the Ha image shows that the 
non-ionizing UV bright stars have migrated out of their 
natal regions. Young, hot stars are more efficient than 
the old stars at heating dust since their emission peaks at 
UV/blue wavelengths where dust is most efficient at ab- 
sorbing photons. Given that both the cold dust and the 
stars most efficient at heating dust are concentrated in 
the spiral arms, it is likely that those stars are the dom- 
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Fig. 6. — The FIR to radio ratio (q) map is shown. This map 
probes the galaxy disk with a physical scale of ^ 645pc. The field- 
of-view of the image is 15' X 23'. 



inant heating source for the 160 /xm dust in the spiral 
arms. 

3.3. Infrared- Radio Correlation 

The far-IR dust emission and the radio continuum 
emission both trace active star formation in galaxies, 
and are known to be tightly correlated o n galactic scales 
l)de Jong et al.llT985t iHelou et al.llT985(l . In Fig. H we 
show the far-IR (42 - 122/xm, FIR) to radio ratio map 
of M81 constructed using the 160 /im resolution MIPS 
images and 20 cm continuum data {q = log(FIR/3.75 x 
10-i2jjz) _ log(S'i.49GHz))- The MIPS bands were com- 
bined to c reate a TIR image and the fractional FIR 
IjHelou et al. 1985) comp onent was estimated per pixel 
using the IPale fc Heloul l)2002() models. The value of q 
ranges from ~1.80 for the AGN dominated nuclear re- 
gion to '-~^2.87 for the bright star forming regions, with 
an intermediate value of ~2.07 for the interarm regions. 
A sim ilar range in q was seen for M33 byfHiDDclcin ct al. 
(|2p03) , although they did not comment on any large 
scale structure seen in q variations as we see in M81. 



Fo r comparison, the average q v alue found for galaxies 
bv lYun. Reddv. fc Gond'cmri|2001[) is 2.34. The variation 
in q is coherent with structures related to the M81 spi- 
ral arms, confirming that the radio-infrared correlation 
is more complex inside individual galaxies than between 
galaxies. This variation will be investigated in detail in 
future papers to explore the d ynamics of cosrnic-ray elec- 
tron propagation and escape IjHelou fc Bicavlll993j) . 

4. CONCLUSIONS 

We presented MIPS images of M81 at 24, 70, and 
160 /xm which reveal a bright nucleus and two well re- 
solved spiral arms studded with bright regions of star 
formation. These images show that M81 has a signif- 
icant amount of cold dust associated with the spiral 
arms. From multiwavelength morphology comparisons, 
the dust heating is argued be dominated by recent star 
formation even at the longest MIPS wavelength. The 
resolved UV and Ha SFRs are always lower than the 
IR SFRs indicating significant dust attenuation, radia- 
tive transfer ef fects, and/ or different stellar ages than 
assumed for the lKennicuttI ^998) SFR calibrations. The 
character of the dust attenuation indicates that the dust 
geometry and/or grain properties are different for re- 
solved regions in M81 than in starburst galaxies. The 
infrared-radio correlation was found to vary by a factor 
of ~6 in the M81 with coherent structures relating to 
the spiral arms. These results illustrates the need for 
additional theoretical and empirical work on how to ac- 
curately combine different SFR indicators, accounting for 
dust, age, and radiative transfer effects, to give an accu- 
rate view of the star formation in regions of galaxies. 
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